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Abstract 

We study the transport property for armchair-edge graphene nanoribbons (AGNRs) with an adatom cou- 
pling to a semi-infinite quantum wire. Using the nonequilibrium Green's function approach with tight- 
binding approximation, we demonstrate that the tunneling current through the system is sensitively depen- 
dent on both the AGNR width and adatom location. Interestingly, when the adatom locates onto a carbon 
atom in the 3 jth chain from the edge of a metallic AGNR, the system shows a transmission gap accom- 
panied by a threshold voltage in I - V curve like a semiconducting AGNR. This effect may be useful in 
scanning tunneling microscopy experimental characterization on graphene samples. 
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Recently, the successful fabrication of graphene in experiment results in intensive attention due 
to its unique properties^ and potential applications in electromechanical resonator,^ p-n junction,"* 
transistor^ and thermoelectric power devices,^ etc. However, due to the absence of gap in the 
energy spectrum, graphene is difficult to be used in the electronic device which requires precise 
control of carrier type and transport behavior. Fortunately, GNRs, finite width of graphene, open 
a gap in the energy spectrum because of the quantum confinement efi'ect,^'^ remedying the draw- 
back of graphene. Thus, GNRs are more suitable promising materials for the design of electronic 
devices. ^"^^ According to the atomic configuration of the edge, GNRs are classified into two types 
of AGNR and zigzag-edge GNR (ZGNR). It has been shown that a pristine AGNR exhibits either 
metallic or semiconducting nature with a gap depending on the number of elementary cells A'^ from 
one edge to another,'^''"* i.e., an AGNR is metallic if N=3p-\ (p is an integer) otherwise {N=3p-2 
or 3/7-3) semiconducting, while ZGNRs always maintain metallic behavior with localized states 
near the Fermi level. 

In this Letter, we study the transport property for AGNRs with an impurity atom absorbed on 
the top of a carbon atom in AGNR. As shown in Fig. 1(a), like the model proposed for STM 
experiment on graphene^^"'^ or surface of materials, a semi-infinite quantum wire (as a probe) 
couples to an AGNR via an adatom. In Fig. 1(b), the (red) dashed line rectangle is the periodi- 
cal unit which is composed of N elementary cells for AGNR,^^'^"* where cell (m,n) contains four 
carbon atoms labeled as fi, y, A and 6 connected by the (blue) thick solid line. Using the NEGF 
approach with the help of tight-binding approximation, we demonstrate that the tunneling cur- 
rent through the system is sensitively dependent on both the width of AGNR and the location of 
the adatom on it. Interestingly, when the adatom locates onto a carbon atom in the 3 7th chain 
from the edge of a metallic AGNR, system shows a transmission gap accompanied by a threshold 
voltage in the current- voltage curve, which is a conspicuous semiconducting transport behavior 
in the semiconducting-AGNR-based system. This effect may be useful in the STM experimental 
characterization on the GNR samples and in the application of graphene-based nanodevices. 

The Hamiltonian for this AGNR-based system reads 

H = Hw + Hg + Hi + Ht, (1) 

where Hw=Tj]Li l^i^](^j~YjJ=i ti(a^jaj+i+h.c.) for the semi-infinite wire with chemical potential pi 
and hopping energy h ; //G=Zm.n,a l^2cl,„,aCm,n,a-'Z<; > t2(c}n,„,aCm',n',a' +h.c.) for the TT-band electron 
of AGNR with chemical potential p2, where <•, •> denotes summing over nearest-neighbor sites 
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with hopping integral t2, Hi=eod^d for the adatom with single energy level eq couples to both site 
7=1 of the wire and lattice site a in AGNR respectively through coupling strength vi and V2 with 
tunneling Hamiltonian Hj={vid^ai + V2d^Cm,n,a + h.c). The potential difference between the wire 
and AGNR (via substrate) is related to the applied bias voltage as [/=//i-/i2- 

Therefore, for the AGNR-based system of an adatom located in cell n, the tunneling current 
can be expressed as 



IniU) 



dc0[fi(C0 - jXi) - f2(,co - /U2)]T„(0J), 



(2) 



where /i(2) is the Fermi-Dirac distribution function for the wire (AGNR), T„ is the n-dependent 

transmission probability which can be obtained by solving full Hamiltonian (1) by means of NEGF 
within the nearest-neighbor tight-binding scheme. 

The isolated retarded Green's function (GF) matrix for the system is defined as 








(3) 





where g'^QQ(a))=(a) - eo + iO^T^ is the GF for the adatom, g[^((o) and g^^„,Q,(<t') are GF for the edge 
site (7=1) of the wire and site a in cell n, respectively. Using the NEGF, one finds^^'^^ 



gii(c^>) = I 



^4tj - (o) - ixi^ 



2t\ 2t\ 

with local density of sates (LDOS) pn{(jj)= -^Im\_g\^{co)] at site 7=1 in the wire and 



(4) 



oj-E^ + iO+ 



(5) 



with LDOS Pm,n,a(oj)=-^Im[g'^„ „ ^t sitc a in the nth cell of AGNR. 

Further, based on the nearest-neighbor tight-binding approximation, one can obtain^^'^"^ the n- 
electron energy spectrum in the low-energy limit as 



3k. 



= U2 ± ?2 A/l - 4cos-^cos— ^ + 4cos2-^, 



(6) 



and the wavefunction 



^m,n,a^kx, kq) — C 



g;<Pofe,*?)sin(^^n), a = j3,y 

eiV»fe,Vsin[A:^(n - 1/2)], a = A,6 



(J) 



under the hard-wall boundary condition with normalized constant C=[2(N + 1)]~ ' , where (pa is 
an arbitrary function and the carbon-carbon bond length has been set to unitary. The wavevector kx 
is continuous and the y-direction wavevector is discretized as kq=qnl{N + \) with subband index 
number q=\,2, - ■ ■ ,N for a perfect AGNR. Hence, the energy gap /S.E between the ^th conduction 
and valence band is given by 2?2|1 - 2cos(A;^/2)|. 

Moreover, the retarded self-energy matrix has four nonzero elements of 2jq(^qj^=vi 
(^0Qf(a0)~^2), which describes the tunneling between the adatom and wire (AGNR), re- 
spectively. Therefore, by employing Dyson equation G'"=[(g'")~^-E'^]~\ one gets the 
full retarded GF for the adatom Gqq{oS) and lastly obtains the electronic transmis- 
sion probability r„=4roi(6o»)|GQQ(6o»)pro2(<^) with linewidth functions roi(a>)=:7r|vippii(6t>) and 

ro2(a>)=7r|y2lVm,n,aM- 

In the following we present some numerical examples for the systems with N=l and 8 AGNR, 
respectively. The potentials fii=-fj.2=U/2 and the hopping energies ti=t2=t, where we have set 
the origin of energy as the Fermi energy {Ep=Q>) at equilibrium case. In our model, when the 
adatom energy eo deviates from Ep along the positive or negative direction, the impurity scattering 
for electronic transmission suppresses the amplitude of tunneling current. But this do not change 
the major characteristics for electronic transmission which is determined mainly by the energy 
spectrum of AGNRs. Thus, to obtain the essence of transport property for our system, we fix 60 at 
Ep and set the coupling strength vi=V2=0.5 without loss generality. 

In Fig. 2, transmission probability r„ versus energy (o (in units of t) at equilibrium case of 
jui(2)=0 and tunneling current /„ (in units of ejh) versus bias voltage U for different adatom posi- 
tions are depicted, respectively. Obviously, the transport property of the system strongly depend 
on the width (or number of cells A^) of AGNR and the adatom location {n) on it. As shown in Fig. 
2(a) for the system with N=l , the semiconducting nature is always kept regardless of n, which is 
signed by a transmission gap A6>~ 0.22 around Ep=Q and a series of signature peaks from sub- 
bands. This gap can be determined from the minimal energy gap AEmin- Consequently, as shown 
in Fig. 2(b), only when the applied bias voltage U is above Aoj/l, carriers in the highest valence 
band (q=5) can jump to Fermi level and take part in transport, thus there exists a threshold volt- 
age Uc^O.ll for the semiconducting system. As U increases, other subbands become propagating 
accordingly one by one, which results in broken-line signs in I-U curve in contrast to the smooth 
I-U characteristic in the previously studied similar graphene-based system. 

In contrast, for A^=8 the system shows a metallic behavior with a zero transmission gap [see 



4 



Fig. 2(c)] and linear I-U characteristics [Fig. 2(d)] when the adatom is located at site yS or 7 in 
cells of n= 1,2,4 because that the sixth (^=6) valence band and conduction band meet at the Fermi 
energy with a zero gap.^^ Interestingly, when the adatom is located at cell n=3, there exists a wide 
transmission gap of A(^«0.54 associated with a large threshold voltage Uc^0.21 for the system 
[see (black) solid lines in Fig. 2(c) and 2(d)]. In this case, the metallic -AGNR-based system 
surprisingly exhibits semiconducting behavior. Actually, this phenomenon can be explained by 
wavefunction (7). As shown by the (black) solid line in Fig. 3(a), for metallic AGNR with N=S, 
the adatom locating at site jS or 7 in cell n=3 results in \i/fm,3j3(y)(kx, h)f'=0, hence the sixth (gapless) 
subband does not contribute to transport in this case. Generally, for any metallic-AGNR-based 
system, the semiconducting transport behavior exists only if the adatom is located at the positions 
of carbons in cells of n = 3 j for a = /3(y) and 3 j - 1 for a = A{6) with integer j, respectively. 
In other words, the zero gap subband is not effective when the adatom is located at the 3 jth 
chain from the edge of AGNR because |^^m,«,aP=0 as ^ = 2(N + l)/3. Comparatively, |(Am,n,ctP 
at the ^'th subband with q'=[2(N - l)/3] (where [• • •] rounds a number to the nearest integer) 
begins to have nonzero value from the Fermi level. Therefore, the threshold voltage Uc is taken as 
t2\l-2cos[q'n/(2N+2)] \ depending on the metallic AGNR width N. Of course, this transition from 
metallic to semiconducting does not exist for either metallic AGNRs transport in the graphene- 
plane^"^^ or graphene-based systems. ^^"^^ 

On the other hand, as an application of this transition effect found here, firstly one can find the 
location of an adatom on an AGNR sample in a STM experiment. Further, for a set of AGNRs 
with same p, the metallic AGNR possesses the widest AE^m when its zero gap subband is not 
absent as shown by the (black) solid line in Fig. 3(b). Therefore, when AEmi,, is small enough to 
be neglected, all AGNRs can be treated as graphene. The turning point is around N=400 which 
corresponds the width of AGNR about 100 nm. Moreover, because of quantum confinement, the 
difference between energy levels in the wider AGNRs is smaller than that of narrower AGNRs. So 
in the wider AGNRs, the number of conduction channels within bias voltage window changes less 
abruptly, which gives rise to a smooth I-U curve as shown with the (black) solid lines in Fig. 4(b) 
and 4(d) despite the existence of weak oscillating in transmission spectrum displayed in Fig. 4(a) 
and 4(c). Therefore, the smooth I-U curve is an visible symbol for the wider AGNRs compared 
to the broken-line signs in I-U curve for the narrower AGNRs [see (blue) dot and (red) dashed 
lines in Fig. 4(b) and 4(d)]. However, a threshold voltage Uc (decreases as increases) always 
exists when n=3 j for a=P(y) and 37 - 1 for a=A(6), respectively, even though is large enough 
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to neglect the finite-size efi'ect [see Fig. 4(d)] in our model. Therefore, one can use the magnitude 
of Uc to identify an AGNR sample in practical STM experiment by our setup shown in Fig.l. 

In conclusion, we have considered an AGNR-based system with a single adatom couples to 
both a semi-infinite wire and an AGNR. Using the NEGF technique within the nearest-neighbor 
tight-binding scheme, it is found that the tunneling current of the system is sensitively dependent 
on both the width of AGNR and the location of adatom on it. Interestingly, when the adatom is on a 
carbon atom in the 3 7th chain from the edge of a metallic AGNR, a metallic-AGNR-based system 
always exists a transmission gap and consequently followed by a threshold voltage Uc in I-U 
curve, which is a conspicuous semiconducting transport behavior in the semiconducting- AGNR- 
based system. Although the values of transmission gap and Uc in real situation may slightly diff'er 
from our theoretical prediction due to the unsaturated edge bonds or other factors which induce 
the variation of subbands for GNRs,^"^^ we believe that the main point here is qualitatively sound. 
And this eff'ect may be useful in the STM experimental characterization on AGNRs and in the 
application of graphene-based nanodevices. 
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List of Figure Captions: 

Fig. 1 : (Color online) (a) Electrons tunneling from a semi-infinite wire to an AGNR on substrate 
via an (red) adatom on a camon site of AGNR. (b) The (red) dashed line rectangle represents the 
mth unit along the x-direction, in which the elementary cell n along the j-direction is composed of 
four carbon atoms (fi, y, X, S) connected by blue line. 

Fig.2: (Color online) Transport property of the system with diff'erent AGNR width and the 
adatom position nifi, y). Left for transmission probability TnipS) and right for tunneling current 
/„(f/), where (a) and (b) for A^=7, (c) and (d) for N=%. 

Fig. 3: (Color online) (a) Squared wavefunctions for different subband indices in metallic 
AGNR of A'^=8 as a function of adatom position, (b) Minimal energy gap as a function of AGNR 
width. Note that the zero energy gap in pristine metallic AGNR {N=3p - 1) is absent in the 
corresponding adatom-attached ribbon as indicated by (black) solid line. 

Fig.4: (Color online) Transport property for metallic-AGNR-based system with different and 
n, where (a) and (b) for n=l, (c) and (d) for n=3. 
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